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Abstract An extension of the supersymmetric U model for correlated electrons is 

given and integrability is established by demonstrating that the model can be 
constructed through the Quantum Inverse Scattering Method using an i?-matrix 
without the difference property Some general symmetry properties of the model are 
discussed and from the Bethe ansatz solution an expression for the energies is 

presented. 



The supersymmetric (SUSY) U model was first introduced in []T[] as an example of a 
system of correlated electrons which is integrable in one dimension as a consequence of 
the Quantum Inverse Scattering Method (QISM) (e.g. see [0]). Such models, which can 
be solved exactly by the Bethe ansatz method, are important in that the exact solutions 
offer non-perturbative results concerning physical behaviour. For the SUSY U model, 
Bethe ansatz solutions have been studied [[3], ||, 0, ||, 0, [J| and several analyses into 
the physical characteristics that the model describes have been undertaken [f|, [5], 0, |T0| ]. 

The construction of the SUSY U model is based on an .R-matrix satisfying the Yang- 
Baxter equation associated with the one-parameter family of minimal typical represen- 
tations of the Lie superalgebra gl(2\l). In terms of the standard notation for electron 
creation, annihilation and occupation operators the local (two-site) Hamiltonian for the 
model reads 

Ki+i = -Y,(^c l+1 „ + h.c.)(i + uyw n *'-° +n ^-^ 

a 

+U{c\ l c\ ] c i+1] c i+n + /i.C.) 
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where U is an arbitrary free parameter. The local Hamiltonian is also invariant with 
respect to the Lie superalgebra gl(2\l) (hence the name). Below, an extension of this 
model will be derived in such a way that integrability is maintained. The local Hamil- 
tonian of the new model reads 

h iji+1 = - E (C 1 - **)cU+i* + h.c.) (1 + [/) 1/2(n - + "' +1 -' ) 
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where now t is an additional free variable which when chosen to be real (along with 
U real and U > — 1) results in a Hermitian Hamiltonian. For the case t = the usual 
SUSY U model is recovered. The extended model bears some similarity with the mul- 
tiparametric SUSY U model constructed in [|11[] but is in fact inherently different. 

The construction of the above model is through the use of a solution of the Yang- 
Baxter equation without difference property in the spectral parameter. It is known that 
the Hubbard model may be derived via the QISM using an i?-matrix which is also with- 
out the difference property [ |12"1 , [13]]. However, for the Hubbard model the Lax operator 



is given as a particular coupling of two auxilliary Lax operators of six- vertex type. The 
construction employed here appears more akin to the generalized chiral Potts mod- 
els given in [[HJ based on representations of quantum algebras at roots of unity For 
these models and the one discussed here the spectral parameters without difference 
property originate from the representation of the underlying algebraic structure. 

In order to demonstate integrability of this model, we begin with the rational limit 
of the U q (gl(2\l)) invariant (i.e. gl(2\l) invariant) solution of the Yang-Baxter equation 
with additional spectral paramters constructed in [ ]X5| , p)| ]. This solution may be written 
in the form 

R{u, (3, a) = — —Pi + P 2 + — -P 3 (2) 

u + a + p u — ol — p — 2 

and satisfies the Yang-Baxter equation 

R 12 (u - v, (3, j)R 13 (u, (3, a)R 23 (v, 7, a) = R 23 {v, 7, a)R 13 (u, p 1 , a)R 12 (u - v, (3, 7). 

Note that this solution of the Yang-Baxter equation has multiple spectral parameters, 
one of which has the difference property while the other two do not. It is worth re- 
marking at this point that the above solution is just one of a plethora of solutions of 
this type which arise naturally as a consequence of the representation theory of the 
type I quantum superalgebras [ |T6i , 

In the above expression for the P-matrix the operators Pj are gl(2\l) invariant pro- 
jection operators. To explain their actions we begin by recalling that gl(2\l) has gener- 
ators Pj, i, j = 1,2,3 satisfying the super commutator relations 

Above, [1] = [2] = 0, [3] = 1. Choosing a four dimensional space with basis {\i) : i = 
1, 2, 3, 4}, a representation ir a of gl(2\l) acting on this space exists with the action of the 
generators given by 
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E\ = \2) (3|, 
Ef = \3) (2|, 
E\ = -\3) (3|-|4) (4|, 
El = -\2) (2|-|4) (4|, 
J E 3 2 = V ^|1)(2| + V / ^TT|3) (4|, 
J E 2 3 = v^|2)(l| + v / ^TT|4) (3|, 
El = -Va |1) (3| + V^TT |2) (4| , 
= -yfa |3) (1| + V^TT |4) (2| , 
£f = a |1> (1| + (a + 1) (|2> (2| + |3> (3|) + (a + 2) |4) (4| . 

Above, the states |1) , |4) are bosonic and |2) , |3) are fermionic. The highest weight 
state is 1 1) with weight (0, 0|a). It is significant here that a is a free complex parameter. 
Considering the tensor product representation 7Tg ® ir a for generic values of a and j3, 
then P 1 projects onto the irreducible submodule with (unnormalized) basis vectors 





= |i>®|i), 








= <Jj3\2) <g> |1) + 




5|2), 


^> 


= y^|3) <g> |1) + 




5|3>, 




= VW + 1)|4>6 


5|1) + A 


/a(a + 1) |1> ® |4)) + V^(|2) ® |3) - |3) ® |2)) 



(3) 

while P 3 projects onto the irreducible space spanned by 

*?) = V^+l) |4) ® |1) + V / /3(/3+l)|l) ® |4)) 
+ y/( a + l)(P+l)(- |2) ® |3) + |3) ® |2)), 

*2) = \/^+^ 2 > ® I 4 ) + v^+T|4) ® |2)), 
= |3) ® |4) + v^TT|4) ® |3)), 

*2> = 14) ® |4> • 

The projector P 2 is obtained from F 2 = I— Pi —P%- From this solution of the Yang-Baxter 
equation we may now construct the transfer matrix 

t{u,(3,a) = trR 0L (u,f3, a)....R 02 (u, f3,a)R 01 (u, f3,a) 



3 



which forms a commuting family in two variables; viz. 

/?,£*), t(v,-y,a)] = 

and thus t(u,(3,a) can be diagonalized independently of both u and In fact, the 
diagonalization of this transfer matrix has already been treated in with the result 
that the eigenvalues are given by 

A(u, /?, a) = ( U ~ a ~i \ A(«, P, a) (4) 



with 
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such that the parameters A, , Uj are solutions of the Bethe ansatz equations 

A^y = g^i, , = ,....,„, 

\Ajfc - a J X k + 1 

fc=l A fc — — 1 j=i — ~ 1 

The i?-matrix possesses the property 

R(0,a,a) = -P 

where P is the Z 2 -graded permutation operator. From here on in we make the param- 
eterization 

(3 = itu + a 

and fix t and a. Writing the i?-matrix now as a function of only u we have that 

R(u = 0) = -P 

and thus by the standard approach of the QISM [g] we can construct a closed periodic 
quantum model where the Hamiltonian is the logarithmic derivative of the transfer 
matrix and expressible as 

L-l 

H = hi(i+i) + h L1 

i=l 



with the two-site Hamitonian given by 

h 



-P. —R(u 

du 



In terms of the projections operators (cf. 
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h = —P 1 ^ ——. 
a a + 1 



u=0 

we have 

-P, - lit P- 



dp 



0=a 



a: 



and an 



which when expressed in terms of Fermi operators leads to ([]]) with U 
overall normalization factor of —2{a + 1) included for convenience. The energies of the 
Bethe states are obtained from the transfer matrix eigenvalues 



E 



s . -idly 
-2(a + 1) A — 
du 
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., „ . a + it Xi 
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In conclusion we discuss some remarkable properties this model possesses. By 
construction, the transfer matrix from which it is derived is gl(2\l) invariant. In fact, 
one may show that the eigenstates obtained by the algebraic Bethe ansatz method are 
highest weight states in complete analogy with other gl(2\l) invariant models studied 
in [ |TB| , |I9L [Z0| ]. However, the local Hamiltonians are not gl(2\l) invariant. Only for the 
global system is the supersymmetry present. Regardless, it is still possible to add ar- 
bitrary chemical potential and magentic field terms to the local Hamiltonian (P which 
do not violate the integrability. 

Construction of the model for L = 2 yields the usual supersymmetric U model, as 
does the construction on an open chain using (an appropriate modification of) Sklyanin's 
approach where the boundary ^-matrices are chosen to be trivial. These unusual 
scenarios are a consequence of the property 



P 



dPo 



dp 



dP 2 



P 



j3=a 



which is precisely the symmetry breaking term for the local Hamiltonians, and the 
presence of which also implies that 



Acknowledgements. 

I wish to thank Angela Foerster, Ruibin Zhang and David De Wit for discussions 
and assistance. This work has been supported by the Australian Research Council. 



5 



References 

[1] Bracken A J, Gould M D, Links J R and Zhang Y-Z 1995 Phys. Rev. Lett. 74 2768 

[2] Faddeev L D, Sklyanin E K and Takhtajan L A 1979 Theor. Math. Phys. 40 194 
Kulish P and Sklyanin E K 1982 Led. Notes Phys. 151 62 

[3] Maasarani Z 1995 /. Phys. A: Math. Gen. 28 1305 

[4] Bediirftig G and Frahm H 1995 /. Phys. A: Math. Gen. 28 4453 

[5] Bariev R, Klumper A and Zittartz J 1995 Europhys. Lett. 32 85 

[6] Hibberd K E, Gould M D and Links J R 1996 Phys. Rev. B 54 8430 

[7] Ramos P and Martins M J 1996 Nucl. Phys. B 479 678 

[8] Pfannmiiller M P and Frahm H 1996 Nucl. Phys. B 479 575 

[9] Pfannmiiller M P and Frahm H 1997 /. Phys. A: Math. Gen. 30 L543 

[10] Jiittner G, Klumper A and Suzuki J 1997 /. Phys. A: Math. Gen. 30 1881 

[11] Foerster A, Links J and Roditi 1 1998 /. Phys. A: Math. Gen. 31 687 

[12] Shastry B S 1988 /. Stat. Phys. 30 57 

[13] Martins M J and Ramos P 1998 Nucl. Phys. B 522 413 

[14] Date E, Jimbo M, Miki K and Miwa T 1991 Commun. Math. Phys. 137 133 

[15] Bracken A J, Delius G W, Gould M D and Zhang Y-Z 1994 /. Phys. A: Math. Gen. 27 
6551 

[16] Delius G W, Gould M D, Links J R and Zhang Y-Z 1995 Int. }. Mod. Phys. A 10 3259 

[17] Delius G W, Gould M D, Links J R and Zhang Y-Z 1995 /. Phys. A: Math. Gen. 28 
6203 

[18] Foerster A and Karowski M 1993 Nucl. Phys. B 396 611 
[19] Links J and Foerster A 1999 /. Phys. A: Math. Gen. 32 147 

[20] Foerster A, Links J and Tonel A 1999 Algebraic properties of an integrable t — J 
model with impurities Preprint |cond-mat/ 99010911 

[21] Sklyanin E K 1988 /. Phys. A: Math. Gen. 21 2375 

6 



